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SECTION  I 


INTRODUCTION 

Among  the  major  advantages  of  laminated  composite  structures  over 
conventional  metal  structures  are  their  comparatively  high  strength  to 
weight  and  stiffness  to  weight  ratios.  As  a  result,  fiber  reinforced 
composite  materials  have  been  gaining  wide  application  in  aircraft  and 
spacecraft  construction .  These  applications  require  joining 
composites  either  to  composites  or  to  metals.  Most  commonly,  joints 
are  formed  using  mechanical  fasteners.  Therefore,  suitable  methods 
must  be  found  to  determine  the  failure  strengths  of  mechanically 
fastened  joints.  A  knowledge  of  the  failure  strength  would  help  in 
selecting  the  appropriate  size  joint  in  a  given  application. 

Owing  to  the  significance  of  the  problem,  several  investigators 
have  developed  analytical  procedures  for  calculating  the  strenqth  of 
bolted  joints  in  composite  materials.  Among  the  recent  studies  are 
those  of  V/aszczak  and  Cruse  I’ll,  Agarwal  r  ?  ] ,  and  Garbo  and  Ogonowski 
['ll.  As  will  be  discussed  in  Section  VII,  the  previous  methods 
provide  conservative  results  and  underestimate  the  failure  strength, 
often  bv  as  much  as  fifty  percent. 

The  major  objective  of  this  investigation  was,  therefore,  to 
develop  a  method  which  a)  predicts  the  failure  strenqth  and  failure 
mode  of  mechanically  fastened  composite  joints  with  better  accuracv 
than  the  existing  analytical  methods  and  b)  can  be  used  readily  in  th» 
design  of  mechanically  fastened  composite  joints.  In  th"  present 
method  first  the  stress  distribution  around  the  hole  is  calculated  by 
the  use  of  a  finite  element  method .  Second,  the  failure  load  and  the 
failure  mode  are  predicted  by  means  of  a  proposed  new  failure 


hypothecs  together  with  Yamada's  [4]  failure  criterion.  On  the  basis 
of  this  analysis  a  computer  code  was  developed  which  can  be  applied  to 
joints  involving  laminates  with  different  ply  orientations,  different 
material  properties,  and  different  configurations,  including  different 
hole  sizes,  hole  positions,  and  joint  thicknesses.  Because  of  the 
accuracy  of  the  method  and  the  flexibility  of  the  computer  code,  the 
code  can  be  applied  to  the  analysis  and  the  design  of  mechanically 
fastened  composite  joints. 


SECTION  II 


PROBLEM  STATEMENT 

Consider  a  plate  ( Lenqth  I,,  width  W,  thickness  H )  made  of  N  fiber 
reinforced  unidirectional  plies.  The  ply  orientation  is  arbitrary, 
but  must  be  symmetric  with  respect  to  the  x^=0  plane  (symmetric 
laminate.  Figure  1).  Perfect  bonding  between  each  ply  is  assumed. 

A  hole  of  diameter  D  is  located  along  the  centerline  of  the  plate 
(x1=0)  at  a  distance  F  from  one  end  of  the  plate.  A  rigid  pin 
(diameter  D),  supported  outside  the  laminate,  is  inserted  into  the  hole 
(Figure  1).  A  uniform  tensile  load  p  is  applied  to  the  plate,  as  shown 
in  Figure  1.  The  load  is  parallel  to  the  plate  (in-plane  loading)  and 
is  symmetric  with  respect  to  the  centerline.  Hence  the  load  cannot 
create  bending  moments  about  either  the  x1 ,  x^  or  x^  axes. 

Moreover,  for  symmetric  laminates,  in  plane  and  bending  effects  are 
uncoupled.  It  is  desired  to  find 

1)  the  stresses  and  strains  in  each  ply, 

2)  the  maximum  (failure)  load  (P  )  that  can  be 

max 

applied  before  the  joint  fails,  and 

3)  the  mode  of  failure. 

Point  2  refers  to  the  fact  that,  according  to  experimental 
evidence,  mechanically  fastened  joints  under  tensile  loads  generally 
fail  in  three  basic  modes  referred  to  as  tension  mode,  shearout  mode, 
and  bearing  mode.  The  type  of  damaqe  resulting  from  each  of  these 
modes  is  illustrated  in  Figure  2.  The  obiective,  listed  in  point  3 
above,  is  to  determine  which  of  these  modes  will  most  be  responsible 
for  the  failure. 

The  calculation  proceeds  in  three  steps.  For  a  aiven  geometry  and 
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SHEAROUT 


BEARING 


2.  Illustration  of  the  three  basic  failure  modes 


load 


1)  the  stress  distribution  around  the  hole  is  calculated, 

2)  the  maximum  (failure)  load  is  predicted  ,  and 

3)  the  mode  of  failure  is  determined. 

The  details  of  these  steps  are  presented  in  Sections  TII--V. 
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SECTION  III 


STRESS  ANALYSIS-GOVERNING  EQUATIONS 

The  stresses  in  the  laminate  are  calculated  on  the  basis  of 
anisotropic  theory  of  elasticity  and  classical  lamination  plate  theory. 
Accordingly,  in  the  analysis  planes  are  taken  to  remain  planes,  the 
strain  across  the  thickness  is  taken  to  be  constant  [e^  ^=f  (  ,  x9  )  1 

and  only  plane  stresses  are  considered  (a  a ^=0 )  .  Under 

these  conditions,  in  the  absence  of  body  forces,  the  condition  of  force 
eouilibrium  can  be  expressed  as  r5] 

9all  +  9a12  =  Q 


3x, 


3x, 


(1) 


9  °21  +  9q22  _  0 


3x. 


3X, 


In  index  notation  en.  (1)  becomes 

a .  .  .  =  0 

1Jr3 


(2) 


o„  is  the  stress  in  the  plane  normal  to  the  x^  axis  and  is  in 
the  Xj  direction.  The  subscripts  i  and  j  may  have  the  values  l  or  ? . 
Consider  now  an  elastic  laminate  of  volume  VQ  containing  a  loaded 
hole  as  shown  in  Figure  3.  Stresses  are  applied  over  the  suface  area 
A- .  The  surface  area  An  is  riqidly  fixed  (no  displacement),  while 
the  surface  area  A  is  free  of  applied  stress.  The  total  surface 

K 

area  is 


*  =  \  +  AR  +  (3) 

Let  us  denote  by  u.  any  arbitrary  displacement 

inside  the  body.  u^  is  a  test  function.  The  only  requirement  is 

that  u.  be  continuous,  differentiable  and  be  zero  on  A^.  By 

multiplying  eq.  (2)  by  iT  and  by  taking  the  volume  integral  of  the 


7 


TOTAL  SURFACE :  A 
LOADED  SURFACE:  AL 
FIXED  SURFACE:  AR 

STRESS  FREE  SURFACE:  AF=A-AL-AR 


Configuration  of  an  elastic  laminate  with 
loaded  hole 


result inq  expression  we  obtain 


///  a • .  .  u. 


dV  =  0 


By  employinq  the  identity 


r  u  =  (a .  .  u .  )  .  —  a..u.  . 

i] ,  j  i  i]  i  i]  x3  i/ 3 


and  by  utilizing  Gauss'  theorem,  eq.  (4)  may  be  written  as 


//  a .  .  n .  u .  dA  -  //  fa .  ■  u  .  dV  =  0 
11  11  1  i  J  J  13  1/3 


where  n.  is  the  unit  vector  normal  to  the  surface.  On  the 
3 

free  surface  A_  the  stresses  are  zero,  while  on  the  surface 

r 

the  displacement  is  zero.  These  conditions  give 


/  [  a  .  .  n  .  u . 
JJ  13  3  i 


dA  =  0 


f  f  a .  .  n .  u. 
n  13  3  i 


.  u .  dA  =  0 


The  force  per  unit  area  (called  surface  traction)  at  each  point 

of  the  surface  area  A^  is  F53 

T .  =  a . .n.  (9) 

i  13  3 

Equations  (6)  -  (9)  yield 

JJ  T±  u.  dA  =  ///  o.j  u.fj  dV  (10) 


The  stress  is  related  to  the  displacement  through  the  stress-strain 

relationship,  which  for  an  elastic  body  is 

a  =  E  e  (ID 

ij  ijkJl  k£ 

The  subscripts  k  and  l  may  take  on  the  values  of  1  or  7. .  The 
strains  are  related  to  the  displacements  u.  by  the  expression 
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(12) 


=  i,!5l  + 

£k  £  2'3x£,  3xk 

By  combining  eqs .  (10)  -  (12)  we  obtain 

///  Eijkt  5i,j  uk,e  dv  -  //  Ti  "i  "  <13) 

Vo  *L 

B .  are  the  moduli  of  elasticity  of  the  laminate.  Because 

of  the  laminate  symmetry,  the  following  simplification  may  be  made 

E.  ..  .  =  E  (14) 

x  jk«.  mn 

The  subscripts  i,j,k,  and  sl  are  related  to  m  and  n  as  follows 


i=j=l  -*■  m=l 
i=  j=2  ->  m=2 
i/j  -y  m=3 


k=£*l  ->  n=l 


k=<.=2  -*•  n=2 


The  reduced  laminate  modulus  E  is  qiven  by 

mm 

N  hp  -p 

P  .  — P 

where  h  is  the  thickness  of  the  p-th  ply.  0^  is  the 

transformed  reduced  stiffness  matrix  for  the  p-th  ply  [6] 

(Appendix  A). 


t 


10 


SECTION  IV 


STRESS  ANALYSIS-FINITE  ELEMENT  METHOD 

In  order  to  perform  the  calculations,  the  problem  shown  in  Fiqure 
1  was  simulated  by  the  geometry  given  in  Figure  4.  Because  of  symmetry 
the  stresses  were  calculated  only  in  one  half  of  the  body.  Along  the 
symmetry  axis,  displacement  is  allowed  only  in  the  direction. 

Along  the  lower  edge  of  the  plate,  displacement  is  allowed  only  in  the 
x^  direction.  The  intersection  of  the  symmetry  axis  and  the  lower 
edqe  is  considered  to  be  rigidly  fixed. 

The  surface  of  the  hole  is  subjected  to  a  surface  traction  T.. 

1 

The  parameter  T^  is  related  to  the  applied  load.  The  spatial 
distribution  of  T^  depends  on  the  magnitude  of  the  applied  lead,  on 
the  material  properties,  and  on  the  geometry  in  a  complex  manner.  It 
is  extremely  difficult,  if  not  impossible,  to  determine  the  exact 
distribution  of  T^  inside  the  hole.  To  overcome  this  difficulty  a 
cosine  normal  load  distribution  was  assumed.  With  this  approximation 
T^  becomes 

T.  =  -  ^  n>  cos  9  (17) 

l  irD  i 

The  anqle  9  is  in  the  x^ -x7  plane  and  is  measured  clockwise 
from  the  x?  axis  (Figure  1).  For  isotropic  materials  the  cosine 
normal  load  distribution  (eg.  17)  was  found  to  represent  closely  the 
actual  load  distribution  C7].  Calculations  performed  by  previous 
investigators  also  showed  that  for  composite  materials  the  stress 
distribution  inside  the  body  is  insensitive  to  the  assumed  load 
distribution  Cl, 3,8].  Therefore,  eg.  (17)  should  suffice  for  the 
purpose  of  the  present  analysis  which  is  to  determine  the  overall 
strength  of  the  joint. 


Equations  (13)  and  (17)  give 

f  f  (  E  . ,  u .  .  u,  .  dV  =  f  f  ~  — ?•  n  •  u  ■ 
'  '  J  ijkJ,  i,j  k,  l  tt  0  l  l 


cos  6  dA 


We  recall  that  are  functions  that  can  be  selected  arbitrarily. 
The  unknowns  in  eq.  (10)  are  the  displacements  u.  .  Once  u. 

li 

are  known  the  stress  at  every  point  can  be  calculated  from 


eqs  (11)  and  (17) 


The  method  of  solution  of  eq.  (10)  is  described 


below. 

Solutions  to  eq.(10)  were  obtained  by  a  finite  element  method 
(FEM).  As  a  first  step  in  the  solution  procedure  the  volume  Vq  is 
subdivided  into  M  subdomains  of  volume  V 

q 

M 

V  =  Z  V  (19) 


Equation  (10)  may  now  be  written  as 


£.  ///  Eijkl  ui,j  uk,l  dV  “  f.  //  -  75  ni  ui  cos6  34 

9=1  Vg  9=1  %  <2°> 


At  is  the  surface  of  an  element  on  the  inside  of  the  hole  where 
Lq 

the  surface  traction  is  applied  ( 0_<  0_< tt/2  )  .  At  any  surface  where  load 

is  not  applied,  A  is  zero.  Advantage  is  taken  now  of  the 

ug 

assumption  that  the  strains  and  e^)  are  independent  of 

the  thickness,  i.e.  the  strains  are  independent  of  and  depend 

only  on  and  x2-  Thus,  the  three  dimensional  grid  consisting  of 

M  volume  elements  may  bo  replaced  by  a  two  dimensional  grid  consisting 

of  M  surface  elements  of  area  s  (Figure  3) 

9 


S  -  L  S 


Equation  (20)  thus  becomes 


Figure  5.  Grid  used  in  the 
hand  figure  is  an 
around  the  hole 


A  l1  Eijk£  Ui,j  uk,£  ds  =  *  I  -  ^  ni  ui  cose  dr 

g_i  g  g_1  rLg  (22) 

F  is  a  line  along  the  boundary  of  the  hole  where  the  surface 

Ij 

traction  is  applied  Figure  5).  Fach  segment  of  this  line 

(denoted  by  r  )  coincides  with  the  boundary  of  an  element  g  adjacent 

to  the  hole.  For  those  elements  which  do  not  lie  along  this  line  r, 

Lg 

eauals  zero.  Isoparametric  4-node  elements  were  used  in  this 

investigation.  The  grid  was  generated  using  a  grid  generator  £9],  The 

grid  sizes  were  unequal.  Smaller  grids  were  used  in  the  vicinity  of 

the  hole  to  obtain  a  better  resolution  of  the  stresses.  Utilizing  the 

symmetry  about  the  x?  axis,  a  grid  (consisting  of  306  elements)  was 

placed  on  one  half  of  the  laminate,  as  illustrated  in  Figure  5. 

The  displacement  in  each  element  can  be  expressed  in  terms  of  the 

displacements  of  the  four  nodal  points  r10l 

u .  =  N  q . 

l  a  ^ia 

(23) 

u .  =  N  q . 

l  a  nia 

The  subscript  a  designates  the  node]  points  (<*  =  1,2,3,  or  4 ) . 

is  the  shape  function  described  in  detail  in  Appendix  B.  q^  is  the 

displacement  at  the  nodal  point  a  in  the  i  direction. 

We  define  a  stiffness  matrix  for  the  g-th  element  as 


i  6ka 


//  E.  ...  N  N  .  ds 
g  ljkS.  a,  l  i,j 

g 


K .9  lg  an  eight  by  eight  matrix.  The  subscript  6  may  take 
13  ka 

on  the  values  1,2,3,  and  4.  The  nodal  d  ispl  acements  ci,  and 

k  a 

g^p  are  independent  of  the  surface  and  line  integrations. 
Accordingly,  egs .  (??)  -  (34)  yield 
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qi1  Kj.3ka  qka  qiS 


M  —  r  4P 

Z  q. 0  /  -  ^  n.  N 

.  ^i3  ttD  i 

1=1  „ 


cos0  dr 


(25) 


Lg 


The  nodal  displacements  q.  are  arbitrary  functions  and  hence 

ip 

eq.  (25)  can  be  written 


v/h 


Ki Ska  qka  Fi3  {26) 

ere  the  global  stiffness  matrix,  K.„.  and  the  load  vector 
3  iska 


F .  are  given  by 

ip 


Ki3ka 


M 
I  K 


g=l 


i3ka 


(27) 


M 


F.  „ 


Z  /  -  ^  n.  N  cose  dr 

g=l  r  *D  1  3 

Lg 


(28) 


The  elements  of  K.  and  the  components  of  the  vector  F.„ 

1  p  KCt  ip 

are  known.  Hence,  g^a  can  be  obtained  from  eq.  (26)  using  the 
Gaussian  elimination  method  [12].  Once  qka  are  known  the 
displacements  u^  are  calculated  from  eq.  (23).  A  computer  code  was 
developed  for  performing  the  calculations  and  for  generating  solutions 
( section  VT ) . 


SECTION  V 


PREDICTION  OF  FAILURE 

In  order  to  determine  the  load  at  which  a  joint  fails  and  the 
mode  of  failure,  the  conditions  for  failure  must  be  established. 

In  this  investigation  the  joint  is  taken  to  have  failed  when  the 
combined  stresses  have  exceeded  a  prescribed  limit  in  any  of  the 
plies  along  an  approximately  chosen  curve  (denoted  as  the  characteristic 
curve) .  The  combined  stress  limit  is  evaluated  using  the  failure 
criterion  proposed  by  Yamada  f4].  The  coordinates  of  the  characteristic 
curve  are  established  by  extending  Whitney  and  Nuismer's  failure 
hypothesis  Cl3]  (developed  for  open,  uni oaded  holes)  to  loaded  holes. 


1 )  Failure  Criterion 

Numerous  criteria  for  failure  have  been  proposed  in  the  past 
[14, 15, 16, 17] .  Although  the  concepts  underlying  the  different  failure 
criteria  may  be  different,  the  results  of  the  various  criteria  are 
generally  quiie  similar.  In  this  investigation  Yamada' s  failure 
criterion  was  adopted  [4].  This  criterion  is  based  on  the  assumption 
that  just  prior  to  failure  of  the  laminate  every  ply  has  failed  due  to 
cracks  along  the  fibers.  The  validity  of  this  assumption  is  supported 
by  tests  performed  during  this  investigation  with  64  ply  graphite-epoxy 
(AS/3501-6)  laminates.  Yamada' s  criterion  states  that  failure  occurs 
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and  are  the  longitudinal  and  shear  stresses  in  a  ply, 

respectively  (x  and  y  being  the  coordinates  parallel  and  normal  to  the 
fibers  in  the  ply) .  X  is  the  longitudinal  tensile  strenqth  of  the  ply. 
Sc  is  the  shear  strength  of  a  symmetric,  cross  ply  laminate  which 
has  the  same  number  of  plies  as  the  laminate  under  consideration.  As 
indicated  in  eq.  (29)  failure  occurs  when  e  is  equal  to  or  greater  than 
unity. 

2 )  Failure  Hypothesis-Characteristic  Curve 

'  ie  hypothesis  is  proposed  here  that  failure  occurs  when  in  any 
one  of  the  plies  the  combined  stresses  satisfy  an  appropriately  chosen 
failure  criterion  at  any  point  on  a  characteristic  curve.  The 
characteristic  curve  (Figure  6)  is  specified  by  the  expression 

rc(e)  =  D/2  +  Fot  +  «oc-Rot)  cos  9  (30) 

The  angle  9,  measured  clockwise  from  the  x?  axis,  may  range 

in  value  from  -tt/2  to  */ 2 .  R  .  and  R  are  the  characteristic 

ot  oc 

lengths  for  tension  and  compression  [13,18].  These  parameters  can  be 
determined  experimentally  by  measuring  the  tensile  and  compressive 
strengths  of  notched  laminates.  RQfc  and  Rqc  depend  only  on  the 
material.  Therefore,  the  coordinates  of  the  characteristic  curve  also 
depend  only  on  the  material,  and  are  independent  of  the  geometry  and 
the  stress  distribution. 

In  this  investigation  the  characteristic  curve  is  used  together 
with  the  Yamada  failure  criterion.  Accordingly  (see  eq.  29),  failure 
occurs  when  the  parameter  e  is  equal  to  or  is  greater  than  unity  at  any 
point  on  the  characteristic  curve 
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No  failure 
Failure 


e  <  1 
ell 


at  r 


r 

c 


(31) 


It  is  emphasized  that  the  above  failure  hypothesis  is  used  here 
in  conjunction  with  the  Yamada  failure  criterion  (eq.  29).  However, 
the  hypothesis  is  general,  and  is  not  restricted  to  Yamada' s 
criterion.  The  characteristic  curve  proposed  here  may  be  used  with 
any  other  failure  criterion. 


3 )  Solution  Procedure 

Whether  or  not  a  joint  fails  under  a  given  condition  is 
determined  as  follows.  For  a  given  load 

a)  the  stresses  ^a\>a2'a\^  are  calculated  in  each  ply 

using  eqs .  (11),  (14)  and  (16),  and  the  FEM  described  in  Sections  IT, 

and  III, 

b)  the  longitudinal  and  shear  stresses  (px»aXy)  are 
evaluated  in  each  ply  employing  the  transformation 


xy 


2  2 

=  a^cos  n  +  c^sin  n  +  2a^2sinTi  cosn 

r^sinn  cosn  +  c^sinri  cosn  +  0^2  (cos^n~sin^n ) 


(32) 


=  -a. 


where  n  is  the  angle  measured  counter  clockwise  from  the  x^-axis 
to  the  x-axis  of  each  ply. 

c)  the  parameter  e  is  calculated  (eg.  29)  along  the 
characteristic  curve 

d)  if  e  equals  or  exceeds  the  value  of  unity  ( e^l )  in  any  ply 
along  the  characteristic  curve,  the  joint  is  taken  to  have  failed. 

The  procedure  outlined  above  is  used  to  predict  whether  or  not 
failure  occurs  under  a  given  load.  Due  to  the  assumption  of  a  cosine 
normal  load  distribution  around  the  hole  (eq.  17),  the  calculated 
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stresses  are  linearly  proportional  to  the  applied  load  P.  This 
fact  together  with  Yamada's  failure  criterion  (eq.  29)  gives 
P  ^  e 

This  relationship  is  utilized  to  determine  the  maximum  load 

(P  )  which  can  be  imposed  on  the  joint.  For  a  given  load  P, 
max 

values  of  e  are  calculated  on  the  characteristic  curve  as  discussed 
above  (points  a-d) .  The  highest  value  of  e  (eQ)  is  then  determined, 
and  the  maximum  load  is  calculated  by  the  expression 


P 


max 


P 

e 

o 


(34) 


The  calculation  procedure  described  in  the  foreqoing  also 
provides  the  location  (annle  0^)  at  which  e  first  reaches  the 
value  of  unity  (e=l)  on  the  characteristic  curve.  (Figure  7).  A 
knowledge  of  9^.  provides  an  estimate  of  the  mode  of  failure.  When 

O  O 

9fis  small  (9f  =  0  )  failure  is  by  the  bearing  mode.  V>Oien  9 ^z45 
failure  is  due  to  shearout.  IThen  0^=c>no  failure  is  caused  by  tension. 
In  summary 


-15°<0  <15° 
30°<ef<60° 
75°  < 0 f  <90° 


bearing  mode 
shearout  mode 
tension  mode 


(35) 


At  intermediate  values  of  9f failure  may  be  caused  by  a 
combination  of  these  modes. 
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curve 


SECTION  VI 


NUMERICAL  SOLUTION 

A  computer  code  (designated  as  BOLT)  was  developed  which  is 
suitabe  for  generating  solutions  to  the  problem  formulated  in  Sections 
II-V.  The  reguired  input  parameters  and  the  output  provided  by  the 
code  are  summarized  in  Table  1. 

A  Fortran  listing  of  the  code  and  a  sample  input-output  are 
included  in  Appendix  C. 
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Table  1.  Input  parameters  required  by  the  computer  code  and  the 
output  provided  by  the  code. 

Input  Parameters 

1)  Ply  properties 

a)  Young's  modulus, 

b)  Shear  modulus,  G ^ 

c)  Poisson's  ratio, 

2)  Ply  orientations 

3 )  Geometry 

a)  hole  diameter,  D 

b)  thickness,  H 

c)  width,  W 

d)  length,  L 

e)  edge  distances,  F 

4)  Characteristic  lengths,  R  ..and  R 

ot  oc 

5)  Longitudinal  tensile  strength  of  each  ply,  X 

6)  Shear  strength  of  cross  ply  laminate  Sc 

Output  Parameters 

1)  Failure  load 

2)  The  failure  mode 

3)  Stresses  in  the  laminate 


SECTION  VII 


RESULTS  AND  DISCUSSIONS 

Results  were  generated  in  order  to  assess  the  validity  and 
accuracy  of  the  method  and  the  computer  code  and  to  compare  the  results 
of  the  present  method  with  the  results  of  other  existing  methods  of 
solutions.  In  addition,  parametric  studies  were  performed  to  evaluate 
the  major  characteristics  of  bolted  joints. 

1 )  I sotropic  and  Orthotropic  Plates 

Stress  distributions  were  calculated  in  isotropic  plates 
containing  both  unloaded  (open)  and  loaded  holes  and  in  orthotropic 
plates  containing  unloaded  holes.  These  problems  were  selected  because 
analytical  solutions  are  available  for  comparisons  with  the  results  of 
the  present  method . 

An  analytical  solution  for  the  stress  distribution  in  an  infinite 
(V7-*-00)  isotropic  plate  containing  an  unloaded  hole  was  given  by 
Timoshenko  [191.  The  stress  distribution  in  such  a  plate  was  also 
calculated  by  the  present  method.  The  parameters  used  in  the  numerical 
calculations  are  given  in  Figure  9.  A  large  width  (w/D=14)  was  used  in 
the  calculation  to  approximate  an  infinite  plate.  The  results  of  the 
present  method  and  the  analytical  solution  of  Timoshenko  are  compared 
in  Figure  8.  There  is  excellent  agreement  between  the  stresses 
calculated  by  the  two  methods. 

The  stresses  in  isotropic  plates  containing  loaded  holes  were  also 
calculated.  Plates  of  infinite  and  finite  widths  were  considered. 
Calculations  were  performed  for  the  parameters  given  in  Figure  9  and 
Table  2.  From  the  calculated  stresses,  the  stress  concentration  factor 
was  determined.  The  stress  concentration  factor  is  defined  as 
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STRESS ,  Oo/cr 


Figure  8.  The  stress  o2  along  the  x-^-axis  in  an  isotropic 
infinite  plate  containing  a  circular  hole. 
Comparison  of  the  present  results  with  the 
theoretical  results  given  by  Timoshenko  [19]. 
Parameters  used  in  the  numerical  calculations: 
o=1.64  MPa,  D=2R=7. 62  mm,  W/D=14,  e/D=14, 

L/D=28 


1 
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The  stress  concentration  factors  obtained  by  the  present  method 


were  compared  to  those  reported  by  previous  investigators  (Table 
2).  The  maximum  difference  in  the  stress  concentration  factors  given  oy 
the  different  methods  is  about  20  percent.  The  stress  concentration 
factor  given  by  the  present  method  differs  from  the  values  given  by 
previous  investigators  at  most  by  15  percent. 

The  stresses  in  an  isotropic  plate  of  finite  width  containing  a 
loaded  hole  are  shown  in  Figure  9.  The  stresses  calculated  by  the 
present  method  are  in  excellent  agreement  with  be  Jong's  approximate 
solution  [8]. 

The  stress  distribution  in  an  orthotropic  plate  of  finite  width 
containing  an  open  (unloaded)  hole  was  also  calculated.  The 
calculations  were  performed  for  a  plate  with  the  symmetric  laminate 
lay  up  of  [0/90]s.  An  analytical  solution  for  this  problem  was 
provided  previously  by  Nuismer  and  Whitney  TlO],  who  modified 
Lekhnitskii's  earlier  solution  r?0"!  for  an  infinite  plate.  The  results 
given  in  Figure  10  show  excellent  agreement  between  the  stresses 
calculated  by  the  present  method  and  by  the  analytical  solution. 

The  aforementioned  comparisons  indicate  that  the  present  method 
predicts  the  stress  distribution  around  loaded  and  unloaded  holes  with 
high  accuracy. 


2' 


Table  2.  Stress  Concentration  Factor  (SCF)  Around  a  Pin  Loaded  Hole 
Contained  in  an  Isotropic  Plate  of  Infinite  Width. 
Comparison  of  Present  Result  with  Those  Obtained  by 
Previous  Investigators. 


Investigators 

SCF 

Present  Result* 

0.985 

Hong  [22] 

0.955 

He  Jong  [8] 

1.058 

Eshwar  et  al  [23] 

0.922 

Bickley  [7] 

r—i 

cc 

• 

c 

♦Present  results  was  calculated  for  the  case: 
D  =  7 .  f>2  mm,  W/D  =  8,  F/D  =  4,  L/O  =  14. 
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STRESS, 


L 


|  Figure  9.  The  stress  along  the  x^-axis  in  an  isotropic 

plate  of  finite  width  containing  a  loaded  hole. 
Comparison  of  the  present  results  with  the 
theoretical  results  given  by  De  Jong  [8]. 

I  Parameters  used  in  the  numerical  calculations: 

D=7 . 62  mm,  W/D=5.0,  E/D=4.0,  L/D=14.0 
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STRESS 


Figure  10.  The  stress  o2  along  the  x^-axis  in  an  orthotropic 
finite  plate  [0/90]  containing  a  circular  hole. 
Comparison  of  the  present  results  with  the 
theoretical  results  obtained  by  Nuismer  and 
Whitney  [18],  Parameters  used  in  the  numerical 
calculations:  Material:  Graphite/Epoxy  T300/5208, 
E1=149.8  GPa,  E2=11.2  GPa,  G12=5.39  GPa,  v12=0.29, 
a=2.3  MPa,  D=24.5  mm,  W/D=3,0,  E/D=4.0,  L/D=14.0 
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2 )  Failure  Strength  and  Failure  Mode 

Failure  strengths  of  mechanically  fastened  composite  joints 
calculated  by  the  present  method  were  compared  to  available  data  and  to 
failure  strengths  predicted  by  other  numerical  methods. 

Failure  strengths  of  joints  in  graphite-epoxy  laminates  were 
measured  by  Van  Siclen  f"211  and  by  Garbo  and  Oganowski  [3  1.  The  test 
conditions  are  summarized  in  Table  3.  The  failure  strengths  of  the 
joints  under  these  conditions  were  also  calculated  by  the  present 
method.  The  material  properties  used  in  the  calculations  are 
summarized  in  Table  4.  Comparisons  between  the  experimental  and 
predicted  failure  strengths  are  given  in  Table  5 .  in  these  tables 
comparisons  between  the  data  and  the  failure  strengths  predicted  by 
Agarwal  f  2  ]  and  by  Garbo  and  Oganowski  r3]  are  also  included. 

As  can  be  seen,  of  the  three  analytical  methods  included  in  Table 
5  the  present  one  predicts  the  failure  most  accurately.  In  most 
cases,  the  failure  load  given  by  the  present  method  agrees  with  the 
data  within  about  10%.  On  the  other  hand,  the  failure  loads  given  by 
the  Agarwal  and  by  the  Garbo  and  Oganowski  methods  are  in  error  by  as 
much  as  50  percent. 

One  point  is  of  interest  here.  When  there  is  a  hiah  fraction  of 
zero  degree  plies  in  the  laminate  (parallel  to  the  body  direction)  the 
laminate  fails  In  the  shearout  mode.  In  this  case  the  calculated 
results  are  sensitive  to  the  value  of  the  laminate  shear  strength 
Sc.  For  example,  when  the  Sc  value  obtained  with  cross  ply 
laminates  is  used  to  calculate  the  failure  load  of  rO/±45/90]s 
laminates  containing  70  nercent  of  0°  plies  the  calculated  and 
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Table  3.  Summary  of  Test  Conditions 
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Table  4.  Material  Properties  Used  in  the  Calculations 


Ply  Properties 

Longitudinal  Modulus  F.^  C>pa(10  ksi) 
Transverse  Modulus  F.  2  Gpa(10  ksi) 
Shear  Modulus  Gpa(10  ksi) 

Poisson  Patio  v-^ 

Tensile  Strength  X  Gpa(10  ksi) 

Shear  Strength  S  Goa (10  ksi) 


T300/SP286 
130  (18.7)r2l] 

8 . 274 ( 1 . 2  )  [21 ] 
5.033(0.73)[71] 

0.30r211 

1 . 23 (0 . 17R  )  [ 21  ] 
0.05(0.0073)r2ll 


AS/3501-6 
130(18.85)0  ) 
13.1(1.9)0] 
5.R6(0.R5)O] 
0. 30[3] 
1.58(0.23)0] 
0.17(0.0l7)r3] 


Laminate  Properties 

Cross-ply  Lamina  :e  Shear 
Strength  Sc  Gpa(10  ksi) 

Characteristic  Length  Rotmm(in) 

Characteristic  Length  Rocmm(in) 


T300/SP2R6 

0.125  (a)(0.0lB) 
1.007 (c)  (d) (e)(0.043) 
3.04R  (e>(0. 12) 


AG/3501-6 

0.204  (b)(0.03) 

0. 12(C)(0.017) 

0.5R4  (0 . 023 ) [ 3  1 

1 .  7 0.068) 


(a)  Tests  with  Glass/Fpoxy  showed  S  to  be  about  2.5  times  higher 
than  the  ply  shear  strength  S  O].  The  2.5  multiplier  was  used 
to  obtain  S  of  Oraphite/Fpoxy  T300/SP2R6  from  the  ply  shear 
strength  S  given  in  [21 7. 

(b)  The  value  of  S  for  AS/3501-6  was  taken  to  be  1.7  times  the 
ply  shear  strength  S. 

(c)  This  value  was  used  for  laminates  containing  more  than  70% 

(by  volume)  of  0  degree  plies. 

(d)  For  T300/RP2R6  the  value  was  chosen  from  [1R]  for  T300/520R. 

(e)  The  values  of  R  for  T300/R°2R6  and  AS/3501-S  were  evaluated 
by  the  method  given  in  [1R]  together  with  the  data  in  [24,25]. 
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Table  5.  Comparisons  Between  the  Experimental  (P)  and  Calculated 
(P  )  Failure  Loads.  Case  lumbers  Correspond  to  Test 
Conditions  Given  in  Table  3. 


Material 

Case  No. 

Lay  Up 

Percent  Di f ference ( 1 -P/Pc ) xlOO 

Present  Results 

Aqarwal ( 1980 ) 

p 

T300/SP286 

1 

ro/±45/90]s 

3.7 

8.7 

• 

2 

ro/±45/9n:s 

0.01 

0.01 

, 

3 

ro/±45/9C]2s 

1.81 

• 

JO 

4 

r0  2/±45/90l s 

0.01 

0.01 

5 

ro/±45/90]s 

8.66 

1.16 

6 

ro/±45/9ms 

11 . 18 

1 . 2 

4 

► 

7 

r±45i2s 

11. 1 

49.43 

8 

ro/9o]2s 

12.3 

49.85 

* 

9 

T02/±45]s 

7.64 

20.0 

I 

Present  Results 

Garbo  and 

>U 

Ogonowski (1981 ) 

A8/3501-6 

10 

ro/±45/9ois 

0.5 

45 

11 

ro/±45/9ois 

3.8 

27 

12 

ro/±45/90ls 

15.0 

14 

4 


4 


4 


4 
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measured  values  differ  by  35  oercent.  On  the  other  hand,  if  the  shear 
stress  of  the  individual  ply  is  used  in  the  calculations  the 
difference  between  the  calculated  failure  load  and  the  data  is  only  15 
percent.  This  indicates  that  the  value  of  Sc  must  be  selected 
carefully. 

Failure  loads  calculated  by  Waszczak  and  Cruse  fl]  are  compared  to 
data  in  Table  6.  Waszczak  and  Cruse's  method  also  yields  failure  loads 
which  are  in  error  by  as  much  as  50  percent.  The  present  method  was 
not  applied  to  Waszczak  and  Cruse's  data  because  the  material 
properties  needed  for  the  calculation  were  unavailable. 

It  is  interesting  to  note  that  the  accuracies  of  all  four  methods 
(present,  Agarwal,  Garbo  and  Ogonowski,  and  Waszczak  and  Cruse)  depend 
on  the  arrangements  of  the  plies  in  the  laminate.  In  general,  the  the 
analytical  predictions  are  most  accurate  for  guasi-isotropic 
laminates,  and  are  least  accurate  for  angle  ply  an^  cross  ply 
laminates.  However,  even  for  angle  ply  and  cross  ply  laminates  the 
present  method  yields  results  within  about  10  percent  accuracy,  in 
contrast  with  the  results  of  other  existing  methods  of  solutions, 
which  may  be  in  error  by  as  much  as  50  percent. 

The  failure  modes  predicted  by  the  present  method  were  also 
compared  to  failure  modes  observed  experimentally.  These  comparisons, 
given  in  Table  7,  show  that  the  present  method  predicts  well  the  mode 
of  failure. 

The  aforementioned  comparisons  between  the  results  of  the  present 
method  and  the  data  show  that  the  method  predicts  with  good  accuracy 
both  the  load  at  which  the  joint  fails  and  the  mode  of  failure.  The 
good  agreements  between  the  predictions  of  the  model  and  the  data 
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Table  7.  Comparisons  of  Predicted  Failure  Modes  with  those  Observed  Kxper imentally 
T-Tension  Mode  R-Rhearout  Mode  R-Rearinq  Mode 
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(illustrated  in  Tables  5-7)  create  further  confidence  in  the  model. 

3 )  Effects  of  Geometry  and  Ply  Orientations 

Parametric  studies  were  performed  to  evaluate  the  effects  of  joint 
geometry  and  ply  orientation  on  the  failure  strength  and  on  the 
failure  mode. 

The  effects  of  joint  width  on  joint  failure  is  illustrated  in 
Figure  11.  In  this,  and  in  subsequent  figures,  the  failure  load  is 
normalized  with  respect  to  the  ultimate  tensile  load  of  the  laminate 
(without  the  hole)  in  the  direction  of  the  applied  load.  As  is  shown 
by  the  results  in  Figure  11,  in  general,  the  maximum  load  the  joint  can 
carry  decreases  as  the  hole  size  decreases,  when  the  width  to  hole 
diameter  ratio  is  greater  than  about  3.  As  the  hole  diameter 
approaches  the  width  (W/D-*!)  the  strength  reduces  to  zero  (P  -*•  0). 

Here  failure  loads  were  not  calculated  for  W /T>  less  then  three  because 
at  such  low  w/b  ratios  the  assumption  of  the  cosine  load  distribution 
(eq.  17)  is  inaccurate  [8]. 

The  effect  of  edge  distance  E  (Figure  1)  on  the  failure  load  is 
shown  in  Fiqure  12.  For  the  lay-ups  [0/± 45/90 !J  and  [0/90] 
increasing  edge  distance  results  in  higher  failure  loads,  as  long  as 
E/D  is  less  than  about  4.  For  higher  edge  ratios  (E/r>>4)  an  increase 
in  edge  distance  does  not  seem  to  influence  siqni f icantly  the  failure 
load.  For  the  lay-up  [027*45  1  ,  the  failure  load  does  not  vary 
significantly  with  E/D.  / 

The  effects  of  ply  orientation  on  the  failure  load  are  given  in 
Figure  13.  This  fiqure  illustrates  the  effects  of  two  parameters  1) 
the  maximum  ply  angle  in  the  laminate  and  7)  the  change  in 
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FAILURE  LOAD, 


T300/SP286 


•  TENSION  MODE 
O  SHEAROUT  MODE 
X  BEARING  MODE 


3  0.3 

Cl 


X - X 


[0/±45/90] 


o - o  — 


[o2/±45] 

Wx  — 

V[0/90]2s 


4  6  8 

WIDTH  RATIO,  W/D 


Figure  11.  The  effects  of  width  ratio  on  the  failure  load  of 
laminates  with  different  ply  orientations.  P^  is 
the  tensile  failure  load  of  laminates  without  holes. 


Parameters  used  in  the  numerical  calculations: 

Material:  Graphite/Epoxy  T300/SP286,  W=38  mm,  E=50.8  mm 
L=203.2  mm,  H=1.067  mm  for  [0/±45/90]s  and 
[0~/±451  ,  and  H=1.18  mm  for  [ 0/9 0  ]  _ 


FAILURE  LC 


FAILURE  LOAD  ( lb f) 
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15  30  45  60 

PLY  CONTINUITY,  A0  (DEGREE) 


Figure  13.  The  effects  of  maximum  ply  angle  <f>  and  ply 

continuity  A6  on  the  failure  load  of  mechanically 
fastened  joints.  Parameters  used  in  the 
numerical  calculations:  Material:  Graphite/Epoxy 
T300/SP286 ,  D=4 . 76  mm,  W/D=5.336,  E/D=2.983, 
L/D=14. 68,  H=l. 397  mm 


FAILURE  LOAD  (N) 


orientation  between  two  adjacent  plies  A9.  The  latter  parameter  is 
referred  to  here  as  "ply  continuity".  The  results  in  Figure  13  show 

4 

that  the  failure  load  increases  both  with  increasing  <f>  and  with 
increasing  A0,  as  long  as  failure  is  by  shearout  mode.  On  the  other 
hand,  the  failure  load  decreases  with  increasing  4>  and  with  increasing 
when  the  failure  is  by  tension  mode.  These  results  indicate  that  care 
must  be  exercised  in  designing  bolted  joints.  If  there  are  no  other 
design  constraints,  the  range  of  ply  orientation  <t>  and  the  ply 
continuity  A0  should  be  determined  with  the  use  of  the  computer  code 
such  that  the  joint  can  withstand  the  highest  load. 


SECTION  VIII 


CONCLUDING  REMARKS 

The  model  and  computer  code  developed  in  this  investigation  can  be 
used  in  the  design  of  mechanically  fastened  joints  involving  fiber 
reinforced  laminates.  The  computer  code  can  be  used  to  determine 

a)  the  optimum  geometry  of  a  joint  for  a  given  load, 

b)  whether  or  not  the  joint  will  fail  under  a  given  load, 

c)  the  failure  load, 

d)  the  mode  of  failure,  and 

e)  the  ply  in  which  failure  first  occurs. 

The  good  accuracy  of  the  method  suggests  that  it  might  be  worth  it 
to  extend  the  method  to  joints  consisting  of  two  or  more  fasteners. 

The  results  of  parametric  studies  performed  with  the  present 
computer  code  show  that  the  material  properties,  joint  geometry,  and 
ply  orientation,  all  effect  significantly  the  strength  of  mechanically 
fastened  joints. 
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Appendix  A  -  The  Transformed  Reduced  Stiffness  Matrix 


P4  PP  22P4 

=  O^cos  n+2  )  sin 'ricos  n+Q22sin  R 

PPP  22P  4  4 

=  (Q11+Q22~4Q66)sin  ncos  n+Q12(sin  n+cos  n) 

P4  PP  22P4 

=  Q^sin  n+2(Q12+2Q6fi)sin  ncos 'n+Q,2cos  n 

PPP  3PPP  3 

=  (Q11-Q12-2Q33)sinncos  n+(Q12~022+2Q33 )sin  ncosn 

ppp  3  PPP  3 

=  ^Q11-Q12-2Q33  ^ sin  ncosn+^0i2-<^22+^®33  ^sinncOS  n 

=  (Q3l+Q22-2Q^*2-2Q33  )sin?ncos2n+Q33(sin4n+cos4n) 
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The  superscript  p  denotes  the  material  properties  of  the  p-th  ply 

p 

and  the  angle  n  is  measured  from  the  x^-axis  to  the  x-axis.  , 

V.y  and  are  the  longitudinal,  transverse  and  shear  moduli  of  the 

P  P 

p-th  ply,  respectively.  v^p  and  are  Poisson's  rations  of  p-th 

ply  and  satisfy  the  relation 


Appendix  B  -  Shape  Function  Used  in  the  Finite  Element  Code. 


In  the  isoparametric  element,  the  geometry  and  the  displacement  of 
the  element  are  described  in  terms  of  the  shape  function  N ffl  by  a 
transformation  froma  Master  Element  in  the  r-s  coordinate  system  to 
the  element  in  the  x-^-x^  coordinate  system  (Figure  14) 


X. 

=  N  (r,s)  x. 

i  = 

1,2 

1 

a  la 

u . 

i 

=  Nc£r's)  ^ia 

a  = 

1,2,3  or  4 

N 

a 

(r,  s)  =  l/4(l+rr#)  (l+ssa) 

-1  < 

r ,  s  <  1 

where  x.  is  the  coordinate  of  node  a  in  the  i-direction  and 
la 

is  the  displacement  of  node  a  in  the  i-direction,  and  ra  and 
the  coordinates  of  node  a  referred  to  th  Master  Element 


’ia 


sa 


are 


Mote  the  property 


12  >"22 


(-1,0) 


X13’X23 


Element  in  x,-x2  coordinates 


Master  element  in  local 
r-s  coordinates 


Figure  14.  Geometry  of  an  element  used  in  the  finite  element 

calculations.  Left:  Element  in  the  x1~x2 

coordinate  system.  Right:  Element  (master 

element)  in  the  local  (r-s)  coordinate  system. 

x.  is  the  coordinate  of  node  a  in  the  i 
la 

direction,  is  the  displacement  of  node  a  in 

the  i  direction  and  (r  , s  )  are  the  coordinates 

a  a 

of  node  a  in  the  r-s  coordinate  system,  i=l,2, 
a=l,2,3  or  4 
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Appendix  C 
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-  Listing  of  the  Computer  Code  "BOLT",  and  a  Sample 
of  Input  and  Output. 
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